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The f(R) gravity models formulated in Einstein conformal frame are equivalent to 
0^ | Einstein gravity together with a minimally coupled scalar field. The scalar field couples 

with the matter sector and the coupling term is given by the conformal factor. We 
apply the holographic principle to such an interacting model in spatially flat and curved 
i | universes. We show that the model leads to a constant ratio of energy densities of dark 

^jy matter to dark energy in a spatially flat universe. In a spatially curved universe, the ratio 

is not a constant and the evolution seems to be model-dependent. However, we argue 
that any cosmologically viable f(R) model can lead to a nearly constant ratio of energy 
densities and therefore alleviate the coincidence problem. 
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It is now generally believed that the universe is currently undergoing a period of accelerated 
expansion. The simplest candidate to produce this cosmic acceleration is the cosmological 
constant, the energy density associated with quantum vacuum. However, there are several 
problems for associating cosmic acceleration with the cosmological constant. First, theoretical 
estimates on its value are many order of magnitude larger than observations [1]. Second, it is 
simply a constant, namely that it is not diluted with expansion of the universe. This latter 
is specifically important in the sense that there are observational evidence [2] demonstrating 
that the cosmic acceleration is a recent phenomena and the universe must have passed through 
a deceleration phase in the early stages of its evolution. This deceleration phase is important 
for successful nucleosynthesis as well as for the structure formation. We therefore need a field 
evolving during expansion of the universe in such a way that its dynamics makes the decelera- 
tion parameter have a signature flip from positive in the early stages of matter dominated era 
to negative in the present stage [3]. There is also another problem which is the focus of the 
present note. It concerns with the coincidence between the observed vacuum energy density 
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and the current matter density. While these two energy components evolve differently as the 
universe expands, their contributions to total energy density of the universe in the present 
epoch are the same order of magnitude. 

As a different point of view, cosmic acceleration may be interpreted as evidence either for 
existence of some exotic matter components or for modification of the gravitational theory. 
In the first route of interpretation one can take a mysterious cosmic fluid with sufficiently 
large and negative pressure, dubbed dark energy. These models are usually invoked a scalar 
field which during its evolution takes negative pressure by rolling down a proper potential. 
In the second route, however, one attributes the accelerating expansion to a modification of 
general relativity. A particular class of models that has recently drawn a significant amount 
of attention is the so-called f(R) gravity models [4]. These models propose a modification 
of Einstein-Hilbert action so that the scalar curvature is replaced by some arbitrary function 



Recently, different models inspired by holographic principle have been proposed to explain the 
cosmic acceleration. The basic idea is that the number of degrees of freedom of a physical 
system scales with its bounding area rather than with its volume [5] . For an effective quantum 
field theory in a box of size L with an ultraviolet (UV) cutoff A, the entropy S scales exten- 
sively as S ~ L 3 A 3 . However the peculiar thermodynamics of black holes has led Bekenstein 
[6] to postulate that the maximum entropy in a box of volume L 3 behaves non-extensively, 
growing as the area of the box. In this sense there is a so-called Bekenstein entropy bound 



Planck mass. It is important that in this relation the length scale L providing an Infrared (IR) 
cutoff is determined by the UV cutoff A and can not be chosen independently. However such 
a non-extensive scaling law seems to provide a breakdown of quantum field theory at large 
scales. To reconcile this breakdown with the success of local quantum field theory in describing 
observed particle phenomenology, Cohen et al. [7] proposed a more restrictive bound. Since 
the maximal energy density in the effective theory is of the order = A 4 , requiring that the 
energy in a given volume not to exceed the energy of a black hole of the same size results in 
the constraint 



If we take the largest value of the length scale L as the IR cutoff saturating the inequality (2), 
we then obtain the holographic dark energy density 



in which 3c 2 is a numerical constant. It is interesting to note that if the length scale L is 
characterized by the size of the universe, the Hubble scale H -1 , then equation (3) gives a 
vacuum energy density of the right order of magnitude consistent with observations [7]. It is 
however pointed out that this yields a wrong equation of state parameter for dark energy, and 
other possible values for L should be chosen such as the size of the future event horizon [8]. 
This conclusion is, however, based on the assumption of an independent evolution of energy 
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densities of dark energy and dark matter. It is shown [9] that, if there is any interaction 
between these two components the identification of L with H" 1 is possible. In particular, the 
authors of [9] [10] argued that such an identification necessarily implies a constant ratio of the 
energy densities of the two components regardless of the details of the interaction. 
In the present note, we investigate the coincidence problem in the context of holographic f(R) 
gravity models. In f(R) models the dynamical variable of the vacuum sector is the metric 
tensor and the corresponding field equations are fourth order. This dynamical variable can be 
replaced by a new pair which consists of a conformally rescaled metric and a scalar partner. 
Moreover, in terms of the new set of variables the field equations are those of General Relativity. 
The original set of variables is commonly called Jordan conformal frame and the transformed 
set whose dynamics is described by Einstein field equations is called Einstein conformal frame. 
The dynamical equivalence of Jordan and Einstein conformal frames does not generally imply 
that they are also physically equivalent. In fact it is shown that some physical systems can be 
differently interpreted in different conformal frames [11] [12]. The physical status of the two 
conformal frames is an open question which we are not going to address here. 
We will work in Einstein conformal frame. The motivation is that in this frame there is a 
coupling between the scalar degree of freedom and matter sector induced by the conformal 
transformation. In this context, we have already studied the coincidence problem without 
any use of holographic principle [13]. It is shown that the requirement of a constant ratio of 
energy densities of the two components, puts some constraints on the functional form of the 
f(R) function. In other terms, there is a class of parameterized f(R) models for which the 
coincidence problem may be alleviated within a particular region of the parameters space. 
Applying the holographic principle to f{R) gravity leads to some important changes in this 
result. We show that identification of IR cutoff with the Hubble scale in a spatially flat universe 
necessarily leads to a stationary ratio of energy densities corresponding to dark energy and 
matter sector regardless of the functional form of the f(R) function. Thus, any holographic 
f(R) model can address the coincidence problem in a spatially flat universe. In a universe 
with a spatial curvature, we will show that the ratio of energy densities is no longer a constant 
and evolves during expansion of the universe. We argue that alleviation of the coincidence 
problem requires conditions that directly link to the cosmo logical viability of f{R) models. 
Let us start with introducing the action for an f(R) gravity theory in the Jordan frame 



where g is the determinant of g^ v and S m is the action of (dark) matter which depends on the 
metric g^ v and some (dark) matter field ip. Stability in matter sector (the Dolgov-Kawasaki 
instability [14]) imposes some conditions on the functional form of f(R) models. These con- 
ditions require that the first and the second derivatives of f(R) function with respect to the 
Ricci scalar R should be positive definite. The positivity of the first derivative ensures that 
the scalar degree of freedom is not tachyonic and positivity of the second derivative tells us 
that graviton is not a ghost. 

It is well-known that f(R) models are equivalent to models in which a scalar field minimally 
couples to gravity with an appropriate potential function. In fact, we may use a new set of 
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variables 

9ixv = & g^u (5) 

0=^lnfi (6) 

where Q = = /'(-R) and /? = This is indeed a conformal transformation which 

transforms the above action in the Jordan frame to the following action in the Einstein frame 
[11] [15] 

Sef = \J dV=$ {jLr - - 2V(<f>)} + S m (9,ue 2WM ^, V>) (7) 

All indices are raised and lowered by g^ u . In the Einstein frame, is a minimally coupled 
scalar field with a self-interacting potential which is given by 

VW)) = ^ (8) 

Note that the conformal transformation induces the coupling of the scalar field <fi with the 
matter sector. The strength of this coupling j3, is fixed to be ^J~^ and is the same for all 
types of matter fields. In the action (7), we take g^" and as two independent field variables 
and variations of the action yield the corresponding dynamical field equations. Variation with 
respect to the metric tensor g^, leads to 

= M- 2 (7* + Tp (9) 

where 

1%, = V M 0V^ - ^,VW 7 - V{4>)g, v (10) 



^ ~ , f-n Sal" ' 1 ' 



-2 5S m (g^ ^) 

are stress-tensors of the scalar field and the matter field system. It is important to note that the 
two stress-tensors T™ and T$ v are not separately conserved. Instead they satisfy the following 
equations 

= -V%t = ^V,0 f m (12) 

We apply the field equations in a homogeneous and isotropic cosmology described by Friedmann- 
Robertson- Walker spacetime 

dv 2 

ds 2 = -dt 2 + a 2 (t){ , - + r 2 {dB 2 + sin 2 6d0 2 )} (13) 
1 — kr z 

where a(t) is the scale factor and k determines the spatial curvature. The parameter k can 
take 0,-1, +1 corresponding to spatially flat, open and curved universes, respectively. We 
take T™ and T$ v as the stress-tensors of a pressureless perfect fluid with energy density p m , 
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and a perfect fluid with energy density p^ = \4> 2 + V(4>) and pressure = \4> 2 — V ((/>). In 
this case, the equations (9) take the form t 

3 b 

3H 2 + — = M; 2 { P<t> + p m ) (14) 

2k 

2H = -M- 2 [(^ + + p m ] + — (15) 

where = ^ is equation of state parameter of the scalar field <f>, and overdot indicates 
differentiation with respect to cosmic time t. We may write the field equations in terms of 
relative densities defined by Q m = ^ and = — where p c = 3M 2 H 2 is the critical density. 
These relative densities then satisfy 

n m + - n k = 1 (16) 

where Qk = -^jp- The conservation equations (12) give 

p rn + 3Hp m = Q (17) 
P<p + 3H{oo <t) + l) P(t> = -Q (18) 

where 

Q = jf^Pm (19) 

is the interaction term. This term vanishes only for = const., which due to (6) it happens 
when /(it!) linearly depends on R. The direction of energy transfer depends on the sign of Q 
or 0. For <fi > 0, the energy transfer is from dark energy to dark matter and for <fi < the 
reverse is true. 

Let us consider time evolution of the ratio r = p m / ' , 

r = 8n - r P± (20) 

P<t> P<t> 

If we combine the latter with the balance equations (17) and (18), we obtain 

f = 3Hr[u+ + (l + ±)^] (21) 

where 

r = — = (22) 

is the decay rate. Now we apply the holographic relation to dark energy density p^ with 
L = H-\ 

P<p = 3c 2 M 2 H 2 (23) 

This gives 

= 6c 2 M 2 HH (24) 



^Hereafter we will use unbarred characters in the Einstein frame. 
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We combine (23) with (15) to obtain 



r + 



^ = -^(l + ^)(r + iy + - (25) 



One can easily check that 

(r + l)c 2 = fi fc + 1 (26) 

which reduces (25) to 

!^( 1 + -^) (Sil + 1) + i 

Substituting this into (24) gives 



H = ~H\l + ^ T )(^ fc + 1) + - (27) 



p = -9c 2 M 2 tf 3 (l + -^ T )(n fc + 1) + 6c 2 M 2 tf A (28) 

When we put the latter together with the holographic relation (23) into the balance equation 
(18), we obtain 

^=3^^-^ < 29 > 

Note that there is no non- interacting limit in our case since T = corresponds to <fi = const., 
or equivalently, ACDM model. We may use (29) in the relation (21) to obtain 

^^^Xji-^r (30) 
For a further step, we introduce the deceleration parameter which is given by 

*=- l -Jp ( 31 ) 
This together with (27) and (29) results in 

From the above expressions, we note that the important quantities u^, r and q which describe 
evolution of the universe appear as functions of T, r and Q k . We will consider two different 

cases : 

1) In a spatially flat universe, Q k = and the relation (30) results in f = or r = constant. 
It is important to note that this result is independent of the details of the decay rate T or the 
interaction Q. Since the interaction is given by the shape of the f(R) function we conclude 
that applying holographic principle to the dark energy density = 3c 2 MpH 2 in a spatially 
flat universe necessarily leads to a constant ratio of energy densities r = p m / P4>, irrespective 
of the form of the f(R) function. The reasoning is simple : from the holographic relation 
(23) one infers that p^ scales like the critical density p c = 3M 2 H 2 . As a consequence, the 
relative density corresponding to must be a constant so that = ^ = c 2 . With this result 
and the fact that in a spatially flat universe Q k = 0, the Friedmann equation (16) results in 
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Q m = 1 — c 2 . Thus p m has the same scaling as p^ and the ratio r is a constant. 
In the flat case, the relations (29) and (32) reduce to 

w « = -( 1 + ;)3^ < 33 ) 
9 = h r - 1/> (34 > 

Contrary to the constancy of r, accelerating expansion requires particular configurations of 
the f(R) function. This is clear from the expression (34) which the requirement that q < 
automatically sets a constraint on the decay rate. 

The relation (34) indicates a signature flip of the deceleration parameter during expansion of 
the universe and evolution of ^. When ^ < r, we have q > and the expansion is decelerating. 
The condition = r defines a transition region in which q = and the universe enters the 
accelerating phase for which > r. It should be noted that the deceleration parameter 
corresponding to a viable f(R) gravity model should exhibit this signature flip in a special 
way so that 

ff>>0 for r = r ,35, 

It means that we should have an accelerating phase after the decelerating one in the expan- 
sion history of the universe. Moreover, the transition region jj — t should lie in recent past. 
These conditions put some constraints on the shape of f(R) functions. This situation is very 
similar to the case presented in [10] when the Hubble scale is taken as the cutoff length. In 
that case accelerated expansion constrains the shape of the interaction term of an interacting 
holographic dark energy model. 

2) Although a flat universe is usually assumed, it is still quite possible that there is a contri- 
bution to the field equations from the spatial curvature [16]. In a spatially curved universe, 
the relation (32) indicates that jj = r remains as the point in which q changes its sign. 
If we combine (30) and (32) and then use the relation (16), we arrive at 

f = 2q(r + 1)# jj^y = 2 1 H ^ (») 

which relates evolution of r to the deceleration parameter q and This implies that f 
changes its sign at the same point that q does, namely at = r. As the relations (30) and 
(36) imply, contrary to the spatially flat case, r is not a constant and evolves during expansion 
of the universe. Although the density parameter Vl^ is still a constant equal to c 2 , due to the 
non-vanishing spatial curvature the Friedmann equation (16) generally gives different scalings 
for p m and p<f,. The rate of change of the ratio r is closely related to the sign of q. In a spatially 
closed universe (Qk > 0) the ratio r is an increasing function of time when the expansion is 
decelerating (q > 0) and it is a decreasing function when the expansion is accelerating (q < 0). 
In a spatially open universe (f2 fc < 0) the reverse is true. 

There are two important contributions to the changes of r, namely, from Q k and q. There 
are some indications [16] that contributions from the spatial curvature should be very small 
(|f2fc| ~ 0.001). On the other hand, changes of r do not also receive strong contributions from 
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q by noting the fact that accelerating expansion of the universe is a recent phenomenon. To 
clarify this point we first note that transition to the accelerating phase takes place when = r. 
Assuming that this condition is satisfied in an appropriately small redshift near the present 
epoch we conclude that we live in an epoch in which The relation (30) then indicates 

that f receives a small contribution from the deceleration parameter. This is consistent with 
the alleviation of the coincidence problem since besides the possibility that the present epoch 
may be a stationary regime at which the ratio r is a constant, it is also quite possible that we 
live in a very special epoch, a transient epoch at which the ratio varies slowly with respect to 
the expansion of the universe *. 

In summary, applying holographic principle to f(R) models in a spatially flat universe neces- 
sarily leads to a stationary energy density ratio independent of the form of the f(R) function. 
In a spatially curved universe, although r does not appear as a constant ratio and details of its 
evolution is model-dependent, all cosmologically viable f(R) models for which the transition 
region jt = r lies in recent past lead to r ps constant and can alleviate the coincidence problem. 
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